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Abstract 

In this paper, preamble-based least squares (LS) channel estimation in OFDM systems of the QAM and offset 
QAM (OQAM) types is considered, in both the frequency and the time domains. The construction of optimal (in 
the mean squared error (MSE) sense) preambles is investigated, for both the cases of full (all tones carrying pilot 
symbols) and sparse (a subset of pilot tones, surrounded by nulls or data) preambles. The two OFDM systems are 
compared for the same transmit power, which, for cyclic prefix (CP) based OFDM/QAM, also includes the power 
spent for CP transmission. OFDM/OQAM, with a sparse preamble consisting of equipowered and equispaced pilots 
embedded in zeros, turns out to perform at least as well as CP-OFDM. Simulations results are presented that verify 
the analysis. 

Index Terms 

Channel estimation, cyclic prefix (CP), discrete Fourier transform (DFT), least squares (LS), mean squared error 
(MSE), orthogonal frequency division multiplexing (OFDM), quadrature amplitude modulation (QAM), offset QAM 
(OQAM), pilots, preamble. 

I. Introduction 

Orthogonal frequency division multiplexing (OFDM) is currently enjoying popularity in both wired and wireless 
communication systems [2], mainly because of its immunity to multipath fading, which allows for a significant 
increase in the transmission rate [23]. Using the cyclic prefix (CP) as a guard interval, OFDM can "reform" a 
frequency selective channel into a set of parallel flat channels with independent noise disturbances. This greatly 
simplifies both the estimation of the channel as well as the recovery of the transmitted data at the receiver. However, 
these advantages come at the cost of an increased sensitivity to frequency offset and Doppler spread. This is due 
to the fact that, although the subcarrier functions are perfectly localized in time, they suffer from spectral leakage 
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in the frequency domain. Moreover, the inclusion of the CP entails a loss in spectral efficiency, which, in practical 
systems, can become as high as 25% [2]. 

An alternative to CP-OFDM, that can mitigate these drawbacks, is provided by a filter bank-based variant 
employing offset quadrature amplitude modulation (OQAM), known as OFDM/OQAM [12]. This scheme builds 
upon a pulse shaping, which is achieved via (a) an IFFT/FFT-based efficient filter bank, and (b) staggered OQAM 
symbols; i.e., real symbols, at twice the symbol rate of OFDM/QAM, are loaded on the subcarriers [21]. This 
allows for the pulses to be well localized in both the time and the frequency domains. As a consequence, the 
system's robustness to frequency offsets and Doppler effects is increased [11] and at the same time an enhanced 
spectral containment, for bandwidth sensitive applications, is offered [1], [22]. Furthermore, although the two OFDM 
schemes can be seen to exhibit similar peak-to-average power ratio (PAPR) performances, the presence of spectral 
leakage in OFDM/QAM may, ultimately, generate higher peak power values [24]. Moreover, the use of a CP is not 
required in the OFDM/OQAM transmission, which may lead to even higher transmission rates [21]Q 

Since nothing is free in this world, the previously mentioned advantages of the OFDM/OQAM come at the cost 
of subcarrier functions being now orthogonal only in the real field, which means that there is always an intrinsic 
imaginary interference among (adjacent) subcarriers [9]. This makes the channel estimation task for OFDM/OQAM 
systems more challenging, compared to OFDM/QAM. OFDM/OQAM channel estimation has been recently studied 
for both preamble-based [13], [15] and scattered pilots-based [9], [14] training schemes. 

The focus of this paper is on the channel estimation task based on a preamble consisting of pilot tones. The 
question of selecting these tones so as to minimize the channel estimation mean squared error (MSE), subject to a 
given training energy, is addressed. The cases of a full preamble, where all subcarriers carry pilots, and a sparse 
preamble, built upon isolated pilot tones embedded in nulls, are separately treated^] It is shown that an optimal 
sparse preamble for OFDM/OQAM can be constructed with L/, equispaced and equipowered pilot tones, where Lh 
denotes the channel length. 

Related results have previously been derived for the case of CP-based OFDM/QAM (CP-OFDM) channel 
estimation. In [19], it is shown that uniform spacing is the best choice given that the pilot tones are equipowered^] 
Equispaced and equipowered pilot tones were shown in [3] to be the optimal CP-OFDM preamble for a given 
training energy that accounts only for the useful signal, excluding the CP. This paper also revisits the problem of 
optimally selecting the pilot tones in CP-OFDM, when the training energy constraint also includes the CP part. It 
is shown that, in this case, the pilots should also be equal. The effects of such a choice on the resulting PAPR are 
also discussed. For optimal CP-OFDM full preambles, it turns out that they can contain simply equipowered (not 
necessarily equal) symbols. A method of constructing such vectors is also developed. In OFDM/OQAM, all equal 
pilots result in optimal full preambles. 

'Nevertheless, this advantage was partly given up in [16] and a CP-based OFDM/OQAM system was proposed for the sake of facilitating 
the data reception process. 

2 These pilot arrangements are also referred to as block-type and comb-type, respectively [5]. 

3 This is no longer valid if there are suppressed (virtual) subcarriers. In such a case, the optimal placement is non-uniform [18]. 
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In the sequel, the task of using extra (more than the minimum required Lh) pilot tones in a sparse preamble 
is considered and it is shown that no extra gain is provided. Furthermore, the case of including data symbols in 
addition to the pilots, in order to save bandwidth, is also considered and it is shown to result in a performance 
degradation for both OFDM systems. Full and sparse preambles are compared and turn out to ultimately yield 
the same estimation performance. The comparison of optimal sparse preambles for CP-OFDM and OFDM/OQAM 
turns out to be generally in favor of the latter. We present simulations results that confirm the theoretical analysis. 

The rest of the paper is organized as follows: In Section [II] we describe the discrete-time baseband equivalent 
model for the OFDM/OQAM and CP-OFDM systems. The way the various preambles are compared, in terms of 
MSE performance, is detailed in Section [Til] Necessary for the paper definitions and results are summarized in 
Section [IV] Section [V] is devoted to the comparative study of the full versus the sparse preamble vectors, for both 
the OFDM/OQAM and CP-OFDM systems. The use of P > Lh pilot tones in a sparse preamble is investigated 
in Section [VI] Our results concerning the various scenarios of including data with the preamble are briefly presented 
in Section I VIII In Section IVIIII the MSE performances of the optimal sparse preambles associated with the two 
systems are compared. An error floor analysis concerning the OFDM/OQAM system is presented in Section [IX] 
Simulations results are reported in Section [X] Section [XI] concludes the paper. 

Notation. Vectors and matrices are denoted by bold lowercase and uppercase letters, respectively. Superscripts T and 
H stand for transposition and conjugate transposition. The complex conjugate of a complex number z is denoted 
by z* . Also, j = \J —\. || ■ || is the Euclidean norm. For a matrix A, (A)ij denotes its entry. The expectation 
and matrix trace operators are denoted by E(-) and tr(-), respectively. I m denotes the mill-order identity matrix, 
while mxn is the all zeros m x n matrix. 

II. System Models 

In this section, basic definitions of the CP-OFDM and OFDM/OQAM system models are presented, along with 
some basic concepts that will be used in the sequel. 

A. CP-OFDM 

Given M subcarriers, the result of the OFDM modulation of a (complex) Mxl vector x is 

s = ^=T H x 

where T is the M x M DFT matrix, with entries (-F)i.j = e - ^ 4 - 7 , i, j = 0, 1, . . . , M — 1. Prior to transmission, 
a CP of length v is prepended to the previous vector, to yield: 



S QAM 



Im 



(1) 



Assume that the CP length is chosen to be the smallest possible one, namely equal to the channel order: v = Lh — 1 
[19]. Moreover, perfect timing and frequency synchronization is assumed. The channel impulse response (CIR), 
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, is assumed to be constant over the duration of an OFDM symbol. The input to 



the OFDM demodulator, after the CP removal, can then be expressed as 

r = ~Hs + w } 

where Ti. is the Toeplitz circulant matrix with the first row given by ho 0% X (m-l ) ^£>>-i ''' ^2 h% 
and w is the noise at the receiver front end and it is assumed to be white Gaussian with zero-mean and variance 
a 2 . The action of the DFT (FFT) then results in 

V = -jfy^r = diag (H , H x ,..., H M _i)x + rj (2) 
v M 



where H m = Yli=o * hie~ 3 ^~ ml , m — 0, 1, . . . , M — 1 is the A/-point channel frequency reponse (CFR) and 
77 = -^=Tw is the frequency domain noise, with the same statistics as w. The CFR estimates, in the least squares 
(LS) sense, can then be computed as 

H m = y^=H m + ^ (3) 

B. OFDM/OQAM 

The baseband discrete-time signal at time instant I, at the output of an OFDM/OQAM synthesis filter bank (SFB) 
is given by [21]: 

M-l 

«OQAm(0 = ^am,nffm,n(0 (4) 

m= n 

where a m .„ are real OQAM symbols, and 

with g being the real symmetric prototype filter impulse response (assumed here of unit energy) of length L g , M 
being the even number of subcarriers, and f m ^ n = <po + f ( m + n ) mod 7r, where ip can be arbitrarily choserQ 
[21]. The filter g is usually designed to have length L g = KM, where K, the overlapping factor, takes on values 
in 1 < K < 5 in practice. The double subscript (•),„,„ denotes the (m,n)-th time-frequency (TF) point. Thus, m 
is the subcarrier index and n the OQAM symbol time index|f| 

The pulse g is designed so that the associated subcarrier functions g m ,n are orthogonal in the real field, that is 




(l)9p,q{l) \ = $m,pSn,q, (5) 



where &y is the Kronecker delta (i.e., 5i.j = 1 if i = j and otherwise). This implies that even in the absence of 
channel distortion and noise, and with perfect time and frequency synchronization, there will be some intercarrier 
(and/or intersymbol) interference at the output of the analysis filter bank (AFB), which is purely imaginary, i.e., 

9m,n( l )9p,q(l) = (6) 



4 For example, in [21], <p m ,n is denned as (m + n)^ — mnir. 

5 The latter should not be confused with the sample time index I. In fact, the temporal distance between two successive symbol instants 
n,n + 1 equals M/2 sample time instants. 
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and it is known as intrinsic interference [9]. Adopting the commonly used assumption that the channel is (approxi- 
mately) frequency fiat at each subcarrier and constant over the duration of the prototype filter [13], which is true for 
practical values of Lh and L g and for well time-localized g's, one can express the AFB output at the pth subcarrier 
and qth OFDM/OQAM symbol as: 

M-l 

Vp.q — Hp qClp q + J S ' S ' H m , n O'm ,n u m'n ~^~Vp,q (^) 
m— n 

(m,n)/(p,g) 



where H pq is the CFR at that TF point, and I PA and r) p<g are the associated interference and noise components, 
respectively. One can easily see that T] p q is also Gaussian with zero mean and variance cr 2 . 

For pulses g that are well localized in both time and frequency, the interference from TF points outside a 
neighborhood Vt PA around (p, q) is negligible. If, moreover, the CFR is almost constant over this neighborhood, 
one can write (0 as 

Up,q ~ Hp^qCp^q + T]p,q (8) 

where 

Cp.q — a p,q + J ^ ' a m,™ w m 9 n (^) 
(rn,ra)Gn Pig 

When pilots are transmitted at (p, q) and at points inside its neighborhood P)g , the quantity in (O can be 
approximately computed. This can then serve as a pseudo-pilot [13] to compute an estimate of the CFR at the 
corresponding TF point, as 

H PA =y™^H p ,q + ^ (10) 

With a well time-frequency localized pulse, contributions to I p q only come from the first-order neighborhood of 
(p,q), namely fi M = {(p± l,g± 1), (p,q± 1), (p± A special case is given by £l p , q = VLp q = {(p± l,q)}. 

This arises when we place three adjacent pilot tones at positions (p ~ 1, q), (p, q), (p + 1, q) and zeros at the rest 
of the first-order neighborhood positions or when we place nonzero pilot tones at all positions in the preamble 
vector and zero vectors around it. If we abuse OQAM modulation (only in the preamble vector) by transmitting 
the complex symbols dp^e 30 , p = 0, 1, . . . , M — 1, 9 € {0, n, ±7r/2}, then, for an arbitrary p, the corresponding 
pseudo-pilot becomes: 

(m,n)enj ig 

which is real. This is because by using the same phase factors e 3 ® in all the subcarriers we get: 

X>m,n(0ffM(l) = (m,n)efi^ (12) 



If the first-order neighbors of (p, q) carry unknown (data) symbols, one cannot approximate the imaginary 
interference in (0. However, by properly choosing one of the neighboring symbols, say at the point (r, s), this 
interference can be forced to zero. Then the pseudo-pilot in (0 becomes real and equal to a p q . The pilot at (r, s) 
is then known as a /zeZp p//of [9]. 
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III. A Fair Comparison Framework 

For the CP-OFDM system, the preamble structure will consist of one complex vector symbol, as it is common 
in the literature [19]. Note that each complex CP-OFDM symbol is equivalent to two real vector symbols in the 
OFDM/OQAM system. We consider an equivalent preamble structure for the OFDM/OQAM system, which consists 
of one nonzero training vector followed by a zero vector symbol. The latter aims at protecting the nonzero part 
of the preamble from the intrinsic interference due to the data section of the frame [13]. Note that in wireless 
standards (e.g., WiMAX [2]), there are sufficiently long guard periods between the uplink and downlink subframes 
and between frames. Thus, there is no need to worry about intrinsic interference on the preamble vector from 
previous frames. Let us make this more clear: 

Definition 1: Let a preamble structure consist of a number of training vector symbols with only one of them 
being nonzero. This nonzero vector will be called the preamble vector. 

Let Ti,T2 be the sampling periods at the SFB outputs for two OFDM systems. Assume that the minimum 
required number of SFB output samples to reconstruct the preamble vector at the receiver is Ri , R2 for each 
system, respectively. Then the following quantity will be needed in making a fair comparison of the two preambles: 

Definition 2: The training power ratio for the preamble structures of two systems is defined as: 

1 gPi 
TPR Pl ' P2 = RlTl 
R2T2 

where £ p i , £ p i are the energies of the corresponding preambles at the SFB outputs in the minimum sample numbers 
i?i,i?2 respectively. 

Note that, if the two systems are of the same type, then T\ = T2 and the sampling periods can be omitted in the 
last definition. We can now clarify what will a fair comparison will be: 

Remark: Assume that p 1; p 2 are different preamble structures in two systems. Then in order to guarantee a fair 
comparison between these two preamble structures, it is necessary that TPR Pl ' P2 = 1. In other words, we will 
require that the systems under comparison spend the same power on the training data at the transmit antenna. If the 
training power ratio is not equal to one, we can scale the output of the SFB for the second preamble by \/TPR Pl ' P2 
to equalize the training powers for the two preamble structures. 

IV. Summary of Definitions and Useful Results 

Throughout this paper, we make the assumption that M/Lh is an integer number, with M being (as usually in 
practice) a power of two. The channel length, L^, will thus also be assumed to be a power of two@ Moreover, 
we will assume (as usual) that we have a sample-spaced channel [6] and that the nonzero part of the CIR is 
concentrated on its first Lh taps. 

Definition 3: By sparse preamble vector we will mean an M x 1 training vector containing Lh isolated pilots 
and zeros at the rest of its entries. 

6 It this is not the case, one may zero pad the CIR to 2^ log 2 L h\ taps. 
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Definition 4: A preamble vector will be said to be full if it contains pilots at all of its entries. 

Definition 5: A preamble vector with L/ 7 isolated pilots and data symbols at the rest of its entries will be called 
sparse-data preamble vector. 

The following results will be useful in the sequel, hence we briefly summarize them here. Their proofs are in 
the appendices. 

Theorem 1: For CP-OFDM, the sparse preamble that minimizes the MSE of the CFR estimates (01, subject to 
a constraint on the energy of both the useful part of the transmitted signal and the CP, consists of equispaced and 
equal pilot tones. 

Proof: See Appendices U HI] ■ 
Remark. Recall that the MSE-optimal sparse preamble for CP-OFDM is built with equispaced and equipowered, 
not necessarily equal pilots, if for the training energy the CP part is not included [19]. However, at least in theory 
this is not fair since the total amount of energy actually spent for training includes the transmission of the CP 
as well. According to Theorem [T] if the preamble optimization is to be performed in a fair way, then the pilots 
should necessarily be all equal. Of course, such a preamble would suffer from a high PAPR. In practice, this can 
be overcome by transmitting equipowered instead of equal tones, at the expense of a small (practically negligible) 
performance loss. 

Although the above result holds true for the corresponding full preamble as well, we also prove that: 
Theorem 2: There are full preambles for the CP-OFDM system that are MSE-optimal subject to a total energy 

budget both on the useful part of the transmitted signal and the CP, which result in equipowered but unequal pilot 

tones. 

Proof: See Appendix HIT1 ■ 
This result indicates that we can construct optimal full preamble vectors that do not suffer from high PAPR values. 
The construction proposed in Appendix HIH is not unique and might be generalized. However, the proposed algorithm 
provides an infinite possible number of such full preamble vectors, verifying that a low PAPR full preamble vector 
construction is possible. 

For the OFDM/OQAM system, the corresponding results are: 

Theorem 3: For OFDM/OQAM, the sparse preamble that minimizes the MSE of the CFR estimates ( [Tol l, subject 
to an energy constraint, is built with equispaced and equipowered pilot tones. 

Proof: See Appendix IIVI ■ 

Theorem 4: Full OFDM/OQAM preambles with all equal pilots are locally MSE-optimal subject to a transmit 
energy constraint. Their global MSE-optimality is assured when the transmit energy constraint is translated at the 
input of the SFB. 

Proof: See Appendix IIVI ■ 
It should be noted that, for the OFDM/OQAM system, the pilot symbols incorporate the corresponding phase factor 
e Wm '". In view of the above results, preamble vectors containing equal symbols are the only or one of the optimal 
solutions. For the sake of analytical convenience, optimal preambles will henceforth be assumed to consist of all 
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equal pilots. 

V. Full vs. Sparse Preamble 

Let Tj Jh = | io + kj^- k = 0, 1, . . . , Lh — 1 j, for a fixed preselected io 6 |o, 1, . . . , — 1 j, be the set of 
indices of the nonzero pilot tones in the sparse preamble. We denote by Tj^ the set {0,1,..., M — 1} \ Tl h of 
the remaining indices. Then, the full preamble vector can be written as: 

x = x L h + X TZ (") 



where x L> is an M X 1 vector with the nonzero pilots at the positions dictated by and zeros elsewhere, 
while x— is the M x 1 vector containing the rest of the pilot tones in the full preamble vector at the positions 
dictated by Xj-^ and zeros elsewhere. It is evident that x^ x— = and ||a;|| 2 = ||a; L J| 2 + \\xj 



—II 2 - 



A. OFDM/OQAM 

The output of the SFB corresponding to the preamble section is then given by (cf. (01 with n = 0): 

M-l M 
soqam(0 = ^ a rn , g m , (l), I = 0, 1, . . . ,L g + — 1 

m=0 

with the nonzero samples located at k = 0, 1, . . . , L„ — 1. Incorporating the phase factors e 3Vm -° into the real-valued 
symbols a m< o, we obtain the complex training symbols x m .o — a m .oe Wm °. Then imposing the restriction that all 
the training symbols are equal, the phase factors can also, without loss of generality, be considered all equal, say 
e 3V . Note that the requirement of equal symbols essentially leads to an abuse of the OQAM modulation, however 
this happens only in the preamble section. This has already been used in [15] in order to enhance the channel 
estimation performance. Using equal symbols of magnitude the last expression is then written as: 



SO qam(0 = W^° QAM 



where g' mn (l) = g (l — n^) e° M ^ 2 ' ■ We may stack the nonzero output samples, corresponding to the 
preamble vector, to get s OQAM = s O qam(0) soqam(I) ••■ soqam (L g - 1) 
where Sqqam = W^ AM \ £ m6lj . <4,o(0) E me x^ <4,o(l) - E me x,„ <L,o " 1) 



T , — 

a Lh 1 „Lh 

~ OQAM ' OQAM' 



and sim- 



ilarly for Sqqam- 

The energy of the SFB output corresponding to the full preamble can be expressed as: 

4qam = HsoqamII 2 - M(l + 2/3)£OQ AM (14) 
while the energy of the sparse preamble as: 

£qqam = II s oqamII 2 = L h £^ AM , (15) 
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as it easily follows from the analysis in Appendix II VI A proper analytical description for the quantity (3, related to 
the intrinsic interference from adjacent subcarriers, will be given later on. The power ratio for the full and sparse 
preamble structures of the OFDM/OQAM system is thus given by 

TPRf , s _ m+m 

lrrt OQAM — r ) 

where the superscript f , s stands for the words /ull over sparse. 

Next, we will evaluate the channel estimation performance of the two preambles in the frequency domain. Let 
us first consider the sparse preamble. The outputs of the AFB at the pilot positions are given by: 

y m ,o = H m . a m . + rim a, m G 1l h 

For the sparse preamble to be meaningful, there should hold that M/Lh > 2. If the pulse is well-localized in 
frequency, we can assume that the received training samples are uncorrelated, since none of them belongs to the 
first-order neighborhood of the rest. Therefore, the maximum likelihood (ML) estimator coincides with the LS 
estimator, and the CFR can be estimated as: 

H m ,o = =H mfi -\ , m£l Lh (16) 

Stacking these estimates in the vector H L and recalling our assumption of a sample-spaced channel, we can obtain 
an estimate of the CIR as: 

h = F7 1 „ r Hj =h + F7 1 vT n'r , 

r i T 

where n' T = n io +M/L h ,o ... ^±ga^M^ , and xr is the L h x L h submatrix of the 

M x M DFT matrix T consisting of its first Lh columns and its rows corresponding to the indices in Tl h - The 

CFR at all M subcarriers can then be recovered by: 



H - F MxL h h H + F MxLh F L l xLh r]' Lh 



where H = 



Ho,o #i,o • • • H 



M _ 1 , and F MxLh denotes the M x Lh submatrix of T, consisting of 
its first Lh columns. Denoting by C Lh the covariance matrix of rj' L , the MSE of the latter estimator is given by: 



MSE OQAM - E 







2" 






H H 







l pi— 1 jjr—H jpH 

MxL h r L h xL h ^L h r L h xL h r MxL h 



By our assumptions, C Lh = £ o% A u I Lh - Additionally, it is known ([19]) that F LhXLh F^ hXLh = F^ hXLh F LhXL 



= Lhl l . Also, using the fact that F^[ xLh F MxLh = MI L , we finally obtain: 

MSE OQAM = OQAM ^ 17 ) 

Special care is needed in the full preamble case. For the assumed construction of the preamble and with a well 
frequency-localized pulse, we can express each output sample of the AFB by: 

Vm.O = H mSi | a m fl + 22 a ™+lfl U m+l,0 I + Vmfi = H„ h0 C„ h o + ?7m,0) TO = 0, 1, . . . , M — 1 

ie{-i.+i} 
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where the same arguments as in the derivation of eq. ( ITTb have been used. In the last equation, for m — the 
value m — 1 corresponds to the (M — l)th subcarrier and for m = M — 1 the value m + 1 corresponds to the 
Oth subcarrier, due to spectrum periodicity in discrete-time. Furthermore, it is easy to show that for any real and 
symmetric prototype function g, we have v%£ +1 = = (3, m = 0, 1, . . . , M — 1, or u™' + 1 = .o = A 

m = 1, 2, . . . , M — 2, and u^'q = —u° ; ^ > _ 1 , u^I^'q = — u o o _1 '°' depending on the value of L g in the factor 
e~iTT m 2 (needed for causality purposes [21]) and due to its dependence on m. Let us consider the first case, 
although handling the second case is equally straightforward. In the light of these facts, we can see that: 



cm.o = \J£ ° QAM (1 + 2/3), m = 0, 1, . . . ,M - 1 
The MSE expression for the full preamble is thus easily seen to be given by: 

MSE ° QAM = £ OQAM^ +2/3)2 (18) 

To make the comparison between the two preambles fair, we have first to equalize the powers at the outputs 
of the SFB's. Scaling the output of the SFB for the sparse preamble by ^TPRqq AM , we achieve this goal. The 
previous analysis holds as is, with the only difference that the MSE for the sparse preamble is now given by: 

MSF S M ° 2 - Lh<j2 (\9) 

C x lr±t OQAM C * V 1 "+" Z P) 

and the ratio of the two MSE's becomes: 

MSE s p QAM = {l+ Lh 
MSE f OQAM M 

Hence, in a dB scale, the sparse preamble is 101og 10 {M/[Lh(\ + 2/3)]} better than the full preamble. 
B. CP-OFDM 

By ©, and due to the complex field orthogonality of the DFT, the ML estimates of the CFR will again coincide 
with the LS estimates and will be as in ( fToT l. The analysis performed for the OFDM/OQAM sparse preamble 



applies also in the CP-OFDM sparse preamble case. Assuming that we transmit pilots of modulus y £$ AM , the 
MSE expression for the sparse preamble will be: 

s _ Ma 2 

M5fi QAM - 7QAM ' 

which coincides with the MSE expression for the full preamble. Once more, to make a fair comparison, we have 
to evaluate the power ratio in this case. Using (Q]) and ([T3l >. the energy transmitted with the full preamble vector is 



£qam — II s qamI| 2 — \\ x \\ 2 + j-7 xH FmxvFmxv x (20) 



CP energy 

where F Mxv is the M x v matrix consisting of the last v columns of T. Clearly, ||a;|| 2 = M£^ AM . To evaluate 
the CP energy, we will need the following: 
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Lemma 1: For the matrix F 



Mxv 



the following relationships hold: 



and 



Proof: For the first relationship, we have: 



For the second one: 



y^X F Mxu F Mxu)i,i — tr (-F Mxv F Mx^j - tr ( F Mxv F MxSj ~ Mv 

i 



/.(^MXf^MxJij — ^-M F Mxv F Mxu^-M — II F Mxi> II 2 — 

where 1 M is the M x 1 all ones vector. The last equation holds because (l/M).F Mxi ,l M represents the last v 
A/-point IDFT coefficients of the rectangular pulse, which are all zero. ■ 
Consequently, since x has been assumed to be of the form xl M , with \x\ = the CP energy in d20l i is 



zero and hence: 



4am = M££ am 



For the sparse preamble: 



£qAM — ll^ijl 2 + ~^7 x L h F Mxi/ F Mxu x L h 

CP energy 



Obviously, 



and for the CP part: 



\\x Li f = L h E^ 



(21) 



(22) 



(23) 



1 1 - 

— ■r H F F H T - _pQAMV^ 

]\,J Lh r Mxi/ r Mxv^L h — ^,j°x 

1=1 



1 " 
= 1 cQAM \ ' 

M ^ 

i=i 



^2 gJ^-m(M-l-i/+0 
e 3T^m{M-l-v+l) 

m = 



= 0, 



(24) 



where we have used the fact that X^m^ 1 e" , ~ m ^ / 1 is zero for any value of M — 1 — v + 1 that is not an 
integer multiple of L^. Since M — + 1 < M — 1 — v + 1 < M — 1, no such multiple exists. 
Thus, the training power ratio for CP-OFDM is: 



TPR f,s _ QAM _ M 
QAM — 75 — Y~ 
QAM 



(25) 



Scaling by <J TPRq^ m the output of the SFB for the sparse preamble, the associated MSE changes to: 

MSF S - ^ 

IVIDJ^qam — f. s „QAM ~ rQAM 

QAM 1 Cx 
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and, finally, the ratio of the MSE's for the two preambles is given by: 



MSE: 



QAM 



MSE QAM 



Lh 

M 



(26) 



The last equation shows that the sparse preamble has a 101og 10 (M/L/J dB better MSE performance than the full 
preamble. However, as it is shown below, this performance difference can be eliminated if the fact that the CIR is 
concentrated on Lh taps is exploited in the estimation procedure. 



C. A Certain Processing to Equalize the Performances 

Consider the previous analysis for the sparse and the full CP-OFDM preamble vectors. The estimates provided 
by the full-preamble vector can be viewed as being of the form: 

H -- 

where e is an M x 1 error vector, of zero mean and covariance I M . Consider the MSE for this preamble: 




mse; 



QAM 



E 



H-H 



.-QAM 



7 QAM 



M 



To convert this estimate to the time domain, we apply to H the transformation I F 



7lH jp 

MxLh MxL f 



MxL h 



If 



we want to bring it back to the frequency domain, we have to apply to the obtained h the transformation F MxLh . 

-i 



This amounts to applying the transformation F MxL 
The MSE is now given by: 

^2 



h ( F MxL h F 



MxL h 



F MxLh to the originally computed H. 



MSE QAM 



QAM 



E 



MxL h 



F MxL h F MxL h 



MxL h < 



7 QAM 



and therefore (|26l ) now simplifies to: 



MSE QAM - MSE QAM 



(27) 



Thus, this kind of processing thus leads to scaling the ratio of the MSE's by M/Lh. In Appendix [Vl we show 
that the corresponding effect for OFDM/OQAM is approximately equivalent to scaling the full preamble MSE by 
Lh(l + 2/3) /M, for practical values of M, Lh- Thus, applying this processing in the full preamble-based estimates 
in the OFDM/OQAM system, we again obtain: 



MSE 



OQAM 



MSE 



OQAM 



Remarks: 



1) F 



MxL h 



(Fm 



MxL h F MxL h 



F MxLh projects H onto the space of M -point CFR's with CIR's of length 
Lh- This has the effect of suppressing the estimation noise in the impulse response tail. In fact, the previous 
processing can be seen to be equivalent with constrained LS [7]. 

2) The previous analysis only holds for sample-spaced channels [20]. 

3) The MSE equivalence of full and sparse preambles is valid only for optimal preambles. 
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VI. Do We Need More Than L h Pilot Tones? 

Let us now consider using a sparse preamble with P > Lh equispaced and equal pilots to estimate the channel and 
compare its performance with that of the sparse preamble containing only Lh pilots as before. For the OFDM/OQAM 
system, this preamble is still a sparse preamble so the optimality of equispaced and equipowered symbols holds 
as well. For the OFDM/QAM system, if the CP length is P — 1, the sparse preamble with equispaced and equal 
pilot tones is the only optimal solution (see App. IIIIIl i. If the CP length is Lh — 1, the sparse preamble with equal 
and equispaced pilot tones is the optimal one as it is verified by the following analysis. The basic assumption is 
that M/P is an integer. Additionally, P < M, since otherwise we end up with a full preamble. Based on our 
assumptions on M, Lh, it is easy to see that P/Lh is an even integer. Also, the placement of the P pilots follows 
the optimal rule, i.e., the nonzero subcarriers belong to any one of the sets {io,io + -p-, ■■■,io + (P ~ 1)^}' 
io = 1, 2, . . . , 4^ — 1, denoted by Ip. 

A. OFDM/OQAM 

One can easily verify (cf. ( [15] )) that the training power ratio for P pilots over Lh pilots is: 

4 OQAM £f=vV(0' 



TPR P ' Ls 
irxl OQAM 



Z^m=0 e 



? OQAM v L s -l of]) v ih_1 P' J "ct m ('" i T J ") 



P 



(28) 



Since M/ P > 2, the spacing of the pilots in the preamble vector is larger than the size of the first-order neighborhood 
in the frequency direction. Relying again on the good frequency localization of the pulses, we can assume that there 
is almost no interference among the pilot symbols. To convert the P CFR estimates to the time domain, we use in 
this case the transformation (^Fp xLh F PxLh ^j Fp xLh , where F PxLh is the P x Lh submatrix of J~ consisting 
of its first Lh columns and its rows corresponding to the indices in Ip. Then, the MSE for the sparse preamble 
with P pilots is given by: 



MSEqqam = tr 



* MxL h y PxL h * PxL h ) * PxLh^P* PxL h [* Px.L h * PxL h j ■* MxL h 



where C P is the analogue of C L in this case. In view of the equal spacing of the P pilots, we have Fp xLh F PxLh = 

2 

PI Lh , and, moreover, Cp — oqam 1 p- Using these results in the last expression, we obtain: 

MSE OQAM = pOQAM ~^ ( 29 ) 

For the sparse preamble with Lh pilots, < fT7b holds and hence, after the power equalization dictated by 

MSEoqam = MSEqqam (30) 

B. CP-OFDM 

For CP-OFDM, the analysis is similar. The MSE's are then given by: 



MSE, 



Ps 

QAM 



Ma 2 L h 
Tqam ~p" 
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and 

MSE F M ^ 



QAM Tpn P,Ls cQAM 
irft QAM C ^ 

For the sparse preamble with P pilots, the energy at the output of the SFB is given by £qam = II II 2 + 
jjXpF mxi>Fmxv x p> where x p is the preamble vector containing P equal pilots at positions dictated by 2p and 
zeros elsewhere. It is easy to see that £q AM = P£^ AM since, with a CP of length v = Lh — 1, 



x pF 'mxvFmxv x p — £^ AM / ] 



=1 



e 3^m{M-l-v+V) 



2 „ p_j 2 



= 4 QAM E 



2=1 



^2 e j¥ m ( M - 1 -' j +') 

m=0 



= 



as before. The training power ratio is easily shown to be TPRq'^ = P/L h in this case as well, and therefore: 

MSEq AM = MSEq AM (31) 

Remark. It thus turns out that using more pilot tones than suggested by the channel length would not result in 
any performance gain. Observe that this MSE equivalence again only holds for optimal (with equispaced and 
equal/equipowered pilot tones) preambles. 

VII. Including Data in the Preamble 

What if the inactive tones in a sparse preamble are employed to carry data symbols? That would help saving part 
of the bandwidth consumed for training. In such a context, and in order to make a fair comparison between the 
preambles, the data power will not be considered as part of the training energy. This is because the data transmission 
is a benefit of the mixed (sparse-data) preamble. With this consideration, the implications of using such a preamble 
in each OFDM system can be easily explored using the previous analysis. We summarize some main results for the 
case when the data symbols have the same modulus with the training symbols. The case of a different pilot-to-data 
power ratio can be similarly handled. 

A. OFDM/OQAM 

For this system, we consider two different preamble constructions, keeping the definition of the sparse preamble 
as we have done so far: 

1 ) Preamble with a nonzero vector followed by a zero side vector: For this preamble, we have to consider two 
cases: 

Scenario 1: The preamble vector contains data at all positions dictated by Ij^. Since the data power is not taken 
into account in the equalization of the powers at the SFB outputs, the sparse-data and the sparse preambles result 
in the same training power at the output of the SFB. However, it can be easily proved for the resulting MSE that: 

MSF sd Ma 2 M 2 V ,„ ,2^ Ma 2 

M&±!> q AM - OQAM + 2^, \ n ™M > pOQAM ™ M;3t OQAM 
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The last formula for MSEqq AM also implies that as a 2 — ► (or SNR — > oo), the sparse data estimate will 
present an error floor, i.e., 

melL h 

Scenario 2: We have data at all positions dictated by except for the positions around the pilot tones. In this 
case, the implicit assumption is that M/Lh > 4. Then, according to our assumptions, the intrinsic interference term 
in the previous scenario disappears and MSEqqam = MSEqqam' 

2) Using a nonzero side vector containing data and help pilots: In this preamble structure, the second column 
contains data at all positions except for the subcarriers in X^ h , which are loaded with the help pilots. Placing the 
help pilots at the same positions with the pilots is justified by PAPR considerations. In fact, it can be easily shown 
that, among the first-order neighbors of a TF point, the strongest intrinsic interference comes from those points 
corresponding to the same frequency. Specifically, 

I m.n I I m,n i i m.n i 
\ U m,n±l\ > l U m±l,nl > l U m±l,n±ll 

Consequently, placing the help pilots at the aforementioned positions leads to help pilots with smaller modulus, 
thus reducing the PAPR. Note that we can place the help pilots at the corresponding positions of the first column 
of the data section. 

The above preamble structure aims at maximally exploiting the OFDM/OQAM system, by placing data at all 
available positions. We can then consider three different scenarios: 

Scenario 1: The pilots are placed at the first column. The side column and the first column contain data at all 
positions except for the positions that belong to the first-order TF neighborhoods of each pilot tone. In this case, 
there is no need to use help pilots. Obviously, the sparse preamble and this preamble use the same training power. 
Since they will both use the same processing to get the channel estimates, they will lead to the same MSE. 
Scenario 2: The first column contains the pilots and data at all other positions except for the positions adjacent to 
the pilots. The side column contains data and the help pilots at the aforementioned positions. Caution is then needed 
in the time durations we need to observe each preamble to collect the training energy. For the sparse preamble, 
this is L g samples, while for the sparse-data preamble it is L g + M/2 samples. After some algebra, it can then be 
shown that the corresponding power ratio is: 

TPRSd , s i^^ OQAM (i + L,(l + 

lrrl OQAM i T r OQAM j , M 

■j—lJh&x 9 ' 2 

where 

E( m+i 
le{±l} \ U m,0 



( m.lV 
{ U m,0 ) 



>0, VtogZl, 



For ( > M/(2L g ), we have TPRqq AM > 1 and hence MSEqqam < mse oqam- Generally, the last inequalities 
hold, especially as L g increases for better TF localization of the prototype function. 
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Scenario 3: This case is similar to Scenario 2, where now we also place data at all positions adjacent to the pilot 
tones in the first column. The power of the help pilots will be even larger than that in Scenario 2, and therefore 

MSEqqam < MSE OQAM in this case as welL 

B. CP-OFDM 

In this context, Xj— is replaced in x by x d , which contains equiprobable, zero mean, constant modulus, 
uncorrelated data symbols at the positions dictated by 2^—. The expected value of the energy transmitted is given 
by: 

E [£qAm] = II^lJI 2 + J~£ x L h F MxyFMxv x L h + Tf-® (^d -^Mxi/^Mxt/^d 

= ^QAM + ~JJ E { x A F Mxv F "lxv X dj (32) 

The term E ^a;|^i r ' Mxi/ i r '^ xly a; d ^ vanishes under a zero mean assumption on the data symbols. Notice that, again 
to be fair, we do not take the power of the data E [||:c d || 2 ] into account. However, the term E (x^ F Mxu-^ P Mxv x dj 
has to be included, since it represents the energy in the CP section due to the data, used by the receiver to 
eliminate the (intercarrier and intersymbol) interference at all, data and pilot, positions. In view of our assumption 
of uncorrelated data, we obtain: 

ij " icn,. h 

= £$ AU (M-L h )(L h -l) (33) 

The power ratio in this case is defined as: 

TPR sd,s _ E [^QAm] _ , , (M - L h )(L h - 1) 

QAM cs - i+ MLl >i W 

QAM N1 

and hence: 

MSE S QAM = < MSEq AM 

QAM 

VIII. OFDM/OQAM Sparse Preamble vs. CP-OFDM Sparse Preamble 

From the previous analysis, it follows that the sparse preamble is generally the best choice for the preamble 
structure, in both OFDM systems. Let us then compare the estimation performances of the two systems, when 
using optimal sparse preambles, and with the same transmitted power for training. Clearly, in both cases the same 
model for the received signal, eq. (O, will hold. Moreover, if the spacing of the pilots, M/Lh, is large enough 
(theoretically equal to or larger than 2), the noise components at the corresponding outputs of the AFB for the 
OFDM/OQAM system will be uncorrelated. If we do not equalize the powers at the SFB outputs, the two MSE's 

£QAM 

will obviously be related as MSEqq AM = oq AM MSEq Am . Defining the power ratio for the two systems, 
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TPRq AM / q AM , as the ratio of the training power in CP-OFDM over that in OFDM/OQAM, and scaling the 
output of the OFDM/OQAM SFB by ^TPR^ AM/OQAM , we end up with: 

£QAM MSEq AM 

1Vi5£j OQAM - c OQAM TPR s 

O x ^ lrrL QAM/OQAM 

We have previously seen that 

cs _ r cOQAM 

c OQAM — ^hO x 

and 

cs _ t cQAM 
£"QAM — ^hO x 

The OFDM/OQAM sparse preamble generates L g nonzero samples at the output of the SFB, while the CP-OFDM 
sparse preamble yields M + v = M + Lh — 1 samples. The sampling rate at the output of the SFB's is the same 
for both systems. Hence, to equalize the energies per time unit for the two schemes, we have to form the power 
ratio as follows: 



rp-p-rjs M+L h — 1' ■• - c, - C2<;\ 

±rrt QAM/OQAM - j_ L - QQAM ~ M + L h - 1 



J h c x 

and finally 



MSEqqam = — ± MSEq AM (36) 



M + Lh-1, 

For example, let Lh — 32. Then, for L g — M, the CP-OFDM sparse preamble turns out to be superior to the 
corresponding OFDM/OQAM sparse preamble, while for L g = KM, with 2 < K < 5, the OFDM/OQAM sparse 
preamble is approximately 3 — 9 dB better. 
Remarks. 

1) Note that (M + Lh — l)/L g is the ratio of the time durations of the transmit pulses employed by the two 
systems. 

2) The performance difference can be even greater if we want to achieve a lower PAPR in the CP-OFDM system. 
We will then have to use unequal equipowered pilots, which leads to a slightly worse performance of the 
CP-OFDM sparse preamble. 

3) Nevertheless, at the cost of increasing the bandwidth in the OFDM/OQAM system, the OFDM/OQAM 
and CP-OFDM sparse preambles can become MSE equivalent, in the following way. Note that, due to 
the good time localization of the OFDM/OQAM pulse, there is always a subinterval of the total pulse 
duration in the OFDM/OQAM system with the same length as the CP-OFDM modulator output, that carries 
almost all of the energy of the pulse. In view of the even symmetry of g, we can consider the subinterval 
[—\(M + L h - l)/2], \(M + L h - l)/2]] around its center, where \a\ denotes the smallest integer that is 
not smaller than a. Then, for practical values of M, Lh, it can be easily verified than: 

\(M+L h -l)/2] 

]T 9 2 (\L g /2] +0«0.99 

!=-[(M+L h -l)/21 
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Therefore, we only need to observe this interval to approximately reconstruct the preamble vector at the 
receiver. Then the transmit pulses in the two systems have approximately the same duration (albeit with 
OQAM bandwidth increased), thus leading to almost the same MSE performance for the two sparse preambles. 
Again, the OFDM/QAM sparse preamble can be slightly better if we use unequal equipowered pilots in the 
CP-OFDM sparse preamble for a lower PAPR. 
4) This last comparison setup does not affect any of the previous results, since we have always compared 
preamble structures for the same system and therefore the same pulse duration. We only need to be careful 
in the last two scenarios of Section I VIII for the OFDM/OQAM system. For Scenario 2, it can be seen that, 
for sufficiently large L g , 

TPT ,sd,s „ (M + L„-l)(l + Q 

li-KoQAM- M + Lh _ 1+ M 

The same result also holds in Scenario 3. This shows that the last comparison setup reduces the MSE 
differences in the sparse-data case. 

IX. Error Floor Analysis for OFDM/OQAM and CP-OFDM systems 

The fact that the intrinsic interference is a part of the error signal dependent on transmit signal components, 
indicates the existence of an error floor behavior of the perfromance curves for the OFDM/OQAM system. In 
CP-OFDM , there is not such a problem due to the orthogonality of the DFT transformation and the use of the 
CP. This orthogonality eliminates the interferences coming from the neighboring symbols on each pilot symbol. 
However, the OFDM/OQAM system possesses orthogonality only in the real field. Therefore, interferences to each 
subcarrier symbol coming from the neighboring subcarrier symbols are inevitable in the presence of a complex CFR. 
Generally speaking, the interference is minimized for large M and small K. Large M leads to better localization in 
the frequency domain, while small K minimizes the number of overlapping OQAM vector symbols in the temporal 
direction. We will analyze the error floor behavior of both systems to prove the aforementioned claims. 

First, note that the process of estimating the channel in N > Lh (N < M) positions in the frequency domain, 
then finding the CIR by translating these estimates to the time domain, and finally obtaining the channel gains at all 
subcarriers through a DFT operation, is essentially a DFT interpolation of the original frequency domain estimates. 
The CFR coefficients, originally estimated through the LS estimator, can be expressed as H N = F NxL h, where 
F NxLh is the N x Lh submatrix of the DFT matrix consisting of its first Lh columns and its N rows corresponding 
to the indices of the frequency domain channel gains we wish to estimate. The final estimates of the frequency 
domain channel gains are therefore given by H = F MxLh (^F^ xLh F NxLh ^j F^ xLh H N . Assuming that N 
is a divisor of M and the N pilots are equispaced, the last expression becomes H — jjF MyiL F NxL H N . 

Consider now the received signal on the mth subcarrier for the OFDM/OQAM preamble: 

y m = H m a m + jH m u m + rj m 

We have dropped the temporal index for notation simplification. Here a m is a real symbol (equal for example to 
± v £° CJAM for a QPSK constellation), u m is an (N — 1) x 1 vector of crosscorrelations of pulses transmitted on the 
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active subcarriers with the pulse transmitted on the rath subcarrier, incorporating the corresponding symbols as well, 
and 7i m an (N— 1) x 1 vector of the channel gains on these subcarriers. There are two ways to estimate the original 
channel gain: Either as H m = J£ = H m +j^^ + ^ or as H m = = iJm _^ +J 2^ + _^_ 

where s m is the sum of the entries of u m coming from the immediately adjacent subcarriers of m. The second way 
to estimate the channel comes from a usual assumption in the OFDM/OQAM system, namely that the channel gain 
can be considered to be constant in the first-order neighborhood of m, especially if M is large and K is small. In 
this case, the received model can be approximated by y m w H m (a m + js m ) + r\ m , which justifies the second way 
of estimating the desired channel gain. In our analysis, and for the case that M/N > 2, we will use the first way 
of estimating the channel, because we do not have to assume anything about first-order neighborhoods or other 
approximations. In this way, our analysis becomes exact. 

Assuming that the noise is zero mean and its variance is a 2 , we obtain \x m = E[H n 



H m +J H ^ Um , where the 



channel is considered to be an unknown but otherwise deterministic quantity and the training symbols deterministic 
quantities. Here E[] denotes the expectation operator w.r.t. the noise statistics. The MSE for the above estimate is 
given by: 



MSE r , 



E\\H„ 



\^m U m\ 



The last equation justifies the existence of an error floor for the OFDM/OQAM system, since, as a 2 



(37) 
-» 0, 



\7i T u i 2 

MSE m — y - — ™>t , i.e., as the SNR increases, the intrinsic interference becomes a dominant phenomenon. 

On the contrary, for the CP-OFDM system, the received signal model is y m = H m a m + rj m and H m = = 
H m + 2m.. Thus the estimate is obviously unbiased and its MSE, given by a 2 jo? m , tends to zero as the SNR 
increases. 

For the case that M/N = 1, i.e., for a full preamble, we will use the second estimate for the OFDM/OQAM 
system. This estimate leads to better performance as it has been shown in [13]. The increase of the magnitude of the 
pseudo-pilot, a m +js m , as opposed to a m compensates for the inaccuracy introduced by considering the channel to 
be constant in every first-order neighborhood, especially in the SNR regime where realistic systems operate. With 

and its 



this estimating method, the mean value of the exact estimate is ji m = E[H m ] = H m - a ^ /s h J— +/ 



MSE MSE m = \H„ 



1 



We may now stack the N estimates of the channel gains for the OQAM system to obtain the N x 1 vector 



H 



H, 



Case 1: M/N > 2 
We then have w 1 — 



nj u i0 n\u %i 

7 2-, 7 1 -. 



% N — 1 'JV-l 



{0, 1, . . . , M/N - 1}. It can be easily shown that: 



and w 2 = 



, where i € 



MSE« = E 



\H-H\\ 



N 2 



P MxL h FNxL h w l 



+ 



N 2 



<7 2 \\F A/ f v r . F 



H 



where w„ 



l 

_a io ' 
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With a 2 — > 0, the error floor results: 

MSE*** = J^\\FMxL h KxL h Wl\ 

Case 2: M/N = 1 

This is the case of the full preamble. Then: 



-jHo 



so 



+ .] 



SM-l 



-jHi 



Sl 



+.]■ 



n T x u x 



, -jHm-i 



and Wn 



where Wn 



od+JSo ' ai+jsi ' ' ' ' ' a M -l+Js M -l 
1 



a&r-i + jsm-1 
Now: 



+ 3- 



ai + jsi ai + jsi 

rp ~] T 



M—l^M—l 
' BM-1 +JSM-X 



MSE 



7lH 



\\ F My.L h F MxL h 



M 2 



o*\\F 



MxL h 1 MxL h 



aM-l+JSM-l 



Thus, the error floor is given in this case by: 
1 



MSE M° r = ^II^MxL^xL^ll| 2 

Remarks: 

1) There is one more source of error floor generation. This is the channel length. If the channel length is too 
large, then the received signal models for the OFDM/OQAM system used in this paper and in the literature do not 
hold any more. We do not assume such a degenerate case in our analysis or in the simulation section. 

2) The above analysis holds only for sample-spaced channels. There is an extra floor generating mechanism if 
the channel is nonsample-spaced. 



X. Simulations 

In this section, we present simulation results to verify our analysis. The channel follows the veh-A model [2]. 
The CIR is initially generated with 29 taps and then zero padded to the closest power of two, that is, Lh = 32 
taps. We plot the normalized MSE (NMSE), i.e., E(\\H - J3"|| 2 /||H|| 2 ), versus the transmit bit SNR (E b /N ). 
The curves are the result of averaging 200 channel realizations. For each channel realization, 300 different noise 
realizations are considered. QPSK modulation is employed. 



A. CP-OFDM 

The results are for M = 1024 subcarriers and a CP length v = 31. Fig. shows the NMSE performance of 
the CP-OFDM system for the full and sparse preambles, where for the full preamble we use the CFR estimates as 
in J3J. We observe that the performance of a sparse preamble with Lh equispaced and equal pilot tones is much 
better. Note that the difference of the performances in theory should be 10 log 10 (M/L^) = 15.05 dB, which can 
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-CP-OFDM Full Preamble 
- CP-OFDM Sparse Preamble 



E h /N (dB) 



(a) 



-CP-OFDM Full Preamble 
- CP-OFDM Sparse Preamble 



E h /N (dB) 



(b) 



Fig. 1. NMSE performance of the CP-OFDM system for the full vs. sparse preamble case: (a) directly measuring the performance in the 
frequency domain; (b) after applying the frequency^time^frequency processing. 



. CP-OFDM 2L h pilots 
_ CP-OFDM L pilots 



E/N(dB) 



(a) 



E/NMB) 



(b) 



. CP-OFDM 2L h pilots 
. CP-OFDM L pilots 



Fig. 2. NMSE performance of the CP-OFDM system for the sparse preamble with 2L^ pilots vs. the sparse preamble with pilots: (a) 
before the power equalization; (b) after the power equalization. 



be seen to agree with the simulation. This difference is independent of the SNR value, which justifies the fact that 
the curves are parallel. Fig.HJ) depicts the result of the processing described in Section [V-CI As expected, the two 
preambles lead then to the same performance. 

Fig. |2] presents the performance of the sparse preamble with Lh pilots versus a sparse preamble with P = 2L/, 
pilots, before and after the power equalization. The latter preamble performs better before the power equalization, 
since it leads to the transmission of more power. After the equalization of the powers, the two preambles perform 
similarly, as expected. 
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Fig. 3. NMSE performance of the CP-OFDM system for the sparse-data vs. the sparse preamble case. 



In Fig. [3] the sparse and sparse-data preambles are compared. We observe that the sparse preamble is better 
when the transmit power is the same in both cases. The difference between the curves can be easily checked to 
approximately follow the theoretical results. For example, for M = 1024 and = 32, the theoretical value of the 
performance difference is 2.87 dB, which shows up in the figure. 

B. OFDM/OQAM 

For the OFDM/OQAM system, and using filter banks given in [17], [4], Fig.[T]translates to Fig. [4] The differences 
between the curves before and after the power equalization can be easily checked to be in accordance with the 
analytical results. For M = 1024, = 32, K — 4 and the adopted pulse g, the theoretical difference before 
the power equalization is 101og 10 {M/[L^(1 + 2/3)]} « 12.5 dB, which can be seen in Fig. [4^. In Fig. [4}}, the 
performances are similar, verifying the result proved in Appendix IVl The performance of the sparse preamble with 
P = 2Lh pilots is compared to that of the sparse preamble with Lh pilots in Fig. [5] For the mixed sparse-data 
case, we choose to implement Scenario 3 as described in Section IVII-A.2I This is the most involved among the 
sparse-data scenarios in the OFDM/OQAM system and an example for this is provided in Fig. [6] Note the error 
floor in the sparse-data scenario. 

C. Comparison 

The sparse preambles for the CP-OFDM and OFDM/OQAM systems are compared in Fig. [7] The superior- 
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M=1024, K=4 
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Fig. 4. NMSE performance of the OFDM/OQAM system for the full vs. sparse preamble case: (a) directly measuring the performance in the 
frequency domain; (b) after applying the frequency— >time^frequency processing. 
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Fig. 5. NMSE performance of the OFDM/OQAM system for the sparse preamble with P = 2L^ pilots vs. pilots, with power equalization. 
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Fig. 6. NMSE performance of the OFDM/OQAM system for the sparse-data vs the sparse preamble, with power equalization. 



ity of the OFDM/OQAM sparse preamble, when the entire transmit pulse duration is considered, is evident. 
The analytical results can be seen to be approximately verified. Thus, for Fig. |7^, the theoretical difference is 
101og 10 \KM/(M + Lh — 1)] ~ 4.5 dB, while for Fig.lTJ), it is approximately 5.9 dB. These values agree with the 
difference of the experimental curves. The two systems, however, perform similarly in the alternative comparison 
setup described in Section [Villi (Remark 3), as shown in Fig. [8] 

XI. Conclusions 

Optimal preamble design for LS channel estimation in CP-OFDM and OFDM/OQAM systems was addressed in 
this paper, for both full and sparse preambles. In contrast to earlier related work on CP-OFDM, the energy spent 
for the CP transmission was also taken into account when assessing the energy budget for training. This turned out 
to lead to the requirement of equal instead of simply equipowered pilot tones for the CP-OFDM sparse preamble. 
Equipowered and equispaced pilot tones were shown to comprise the optimal sparse preamble for OFDM/OQAM. 
Possible gains from loading data on the inactive subcarriers of a sparse preamble were also investigated. The sparse 
preamble with as many pilot tones as channel taps turned out to be generally the best choice in terms of both 
estimation performance and economy. The OFDM/OQAM optimal sparse preamble was compared with that of 
CP-OFDM and shown to allow for a significantly better performance, provided the whole pulse is transmitted when 
training. Nevertheless, it will perform similarly to CP-OFDM, at the cost of bandwidth expansion, if the tails of the 
(well time-localized) pulse are left out in the transmission of the preamble. Apart from the bandwidth difference 
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Fig. 7. NMSE performance of the CP-OFDM and OFDM/OQAM sparse preambles: (a) M = 512, K = 3; (b) M = 1024, if = 4. 
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Fig. 8. NMSE performance of the CP-OFDM and OFDM/OQAM sparse preambles with the alternative comparison setup: (a) M = 512, 
K = 2; (b) M = 1024, K = 4. 



between the two systems, this shows the fundamental structural similarity of the two systems. Our analytical results 
were confirmed via simulations. 

Appendix I 

CP-OFDM: Sparse Preamble With Equispaced and Equal Pilot Tones 

It was proved in [19], [3] that a sparse preamble of Lh pilots is MSE-optimal subject to a training energy 
constraint when it is built with equipowered and equispaced pilot tones. That energy constraint did not include the 
energy spent for CP. Our goal here is to determine the optimal sparse preamble when the training energy includes 
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the CP energy as well. 

Let I[ Jh = {io, ii, . . . , iz, h -i} be the set of indices of the nonzero pilot tones in the sparse preamble, and denote 
by T-£- the set {0, 1, . . . , M — 1} \lh h , consisting of the indices of the null tones. Stacking the Lh CFR estimates 
in the vector H T , we can find the CIR as h = FJ 1 ^ , H T = h + -F7 1 v r s, where e is the Lh x 1 vector with 
entries rj m /x m , m G Ih h , and F L xL is the Lh x Lh submatrix of the M x M DFT matrix T consisting of its 
Lh first columns and its rows corresponding to the indices in T^ h . Thus, the MSE of the above estimate is given 
by 



MSE Lh = tr 



H 



^L h \ r L h xL h r L h xL h 



(38) 



where C Lh = er 2 diag ^l/|xi | 2 , l/|xij 2 , . . . , l/\xi Lh _ 1 | 2 ^j is the covariance of e. Without loss of generality, we 

2 2 2 

can order the diagonal elements of C r in descending order, as follows , - p > , CT 1 2 > • • ■ > 1 — - — nr, and the 
eigenvalues of (^F LhXLh F^ hXLh ^j in ascending order, i.e., Ao < Ai < • ■ ■ < Then [?, Lemma 1]: 

L h -1 n 

- E uTT2 X ™ ( 39 ) 



MSEr . = tr 



L h xL h r L h xL 



L 

H 



m=0 



with equality if and only if the matrix in the brackets is diagonal. Let the eigenvalue decomposition (EVD) of 
F LhxLh^LhxLh^j be UAU H , with U being its LhxLh eigenvector matrix and A = diag(Ao, Ai, . . . , Xz h -i). 
To optimize the selection of F LhXLh , we must find the optimal placement of the pilots in the training vector. This 
choice will determine U and A. In order to satisfy ( 139b with equality, we have to make a placement that yields 
U = I Lh . It is known [19], that a placement of the pilot tones with this property is the equidistant one. We will 
focus on this placement now to find its optimal pilot tones. In the next appendix, we show that the sparse preamble 
just obtained is the globally optimal one. 

After the CP insertion, the energy reaching the transmit antenna is given by: 

I! S QAM II 2 = IMP + ~j^ xH F MxuFmxv x 

Lh — 1 -, L h — 1 v 



/ j 1 m x ik Xi ' 



E^ 1 



m— fc,m— l — l 

Thus, the optimization problem we have to solve is stated as: 

a 2 Lh ~ 1 1 

min — V pr (40) 

such that (s.t) ^ \x lm \ 2 + jj ^ x* k x im e J ^ {m - k){M - 1 - /+l) < £ (41) 

m=0 k,m=0 1=1 

where 2"£ ( is any (yet fixed) of the M/Lh equispaced placements of the pilot tones and £ the total transmit energy 
available for training. 

Proposition 1: An optimal solution for the problem d40t . (l4lT i is given by equal symbols. This yields a local 
minimum of the constrained problem. 
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Proof: Forming the Lagrangian function for the above problem, we get: 



L h -1 



a 1 

J (xi , Xi ± , . . • , Xi Lh _ 1 ) = — y ] j— 

rB: 



L h -1 

El 

m=Q 



^ E <»<™E 

k.m— 



7-^(m-fc)(M-l-„+0 



A 4 



= 0, m=Q,l,...,Lh-l 

^h-i i 



where is the Lagrange multiplier|_| Setting the gradient of J to zero, we obtain: 

^ + ^E 1 <E^ (m " )(M - 1 - !/+0 

fc=0 i=l 

Multiplying by Xi m and summing over to, /i can be found as fi = X^m'S) 1 i^ 1 pr , where we have used the 
constraint with equality. 

Consider the preamble vector with equal training symbols. Then by the energy constraint we can easily find that 
|xi m | 2 = \x\ 2 = jj-, since the CP energy part vanishes (cf. d24]>). We can check that the Lagrange equations are 
jointly satisfied by this point: 



dj 



dxi 



o" . / 1 

tX 



1 



£ Lh\x\ 2 







(x,x,...,x) "h\V\- C \X\' \X\' \C U h \. 

Furthermore, it can be easily checked that the Hessian in a neighborhood of (%,%,..., x) is always positive definite. 
Hence, the equal training symbols lead to a local minimum in our optimization problem. 
To show this, we have: 



d 2 j 



dxi dx* 



Lho* 



MS 2 

(x ,x .. . .,x) l—l 



E^ 



]ft(m-q)({M-l)-v+l) 



2£ 2 



and 



d 2 J 



dx im dx* 



3ML h + 2L\ - 2L h - M a 2 
M 2E 2 



(a:,a:,...,a:) 

For an arbitrary vector y in the feasibility set of our optimization problem, we can write: 

* 9 2 J 

2^fc,m=0 VkVm r, * \(x,x,...,x) 



dx ik dx* rn 

■> Lh-1 



MS 2 

2£ 2 

Lh<7 

MS 

f_ 

2£< 



E E ' 





1=1 






Lh-l 

E y k Vm ' 


+ 2L\ - 


2L h 


- M a 2 


f M 




2£ 2 


2 L h -1 v 

2 E ffc^E 

k,m=0 1=1 


e j^(m-fc)((Af-l)- 


-v+l) 


Lh{Lh 
M 


Lh-1 

E ^™ ■ 

/c,m— 0,fc^m 


3ML h + 2L\ - 


2L h 


— M a 2 


\ - — 

M 




2£ 2 



Lk-1 



\Vk\ 



k=0 



Lk-1 



E iw* 



Lh-1 

E 

fe=0 
2 



\Vk\ 



k=0 



7 In the complex field, we should consider the real part of the product of the Lagrange multiplier and the constraint. Nevertheless, in this 
problem, due to the symmetry of the constraint set, it can be proved that both approaches lead to the same result. 
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Furthermore, 



L h -l 

E VlVn 




= (L h -1) \Vk\ 2 <(L h -l)£ 



k=0 



where in (a) we have used the Cauchy-Schwartz inequality, thus 



a 
'28* 



L h -1 

E VlVm > 



(L h - l)a 2 
2£ 2 



E 

fe=0 



\Vk\ 



We can write: 



L h -1 



M£ 2 



E ^»E 



r 2 " 

o j|^(m-fc)((M-l)-r/+i) _ \ - 

M£ 2 ^ 



fe,m=0 1=1 

The positivity of the Hessian is satisfied if 



L h -1 



E N 



E a 

?ri=0 



it- «■ 



> 



MIL h + 2L\ - 2L h - M L h (L h - 1) L h - 1 \ a 



2M 



M 



£2 -L h ^>0 



which holds for any Lh > 1. 



The resulting (time-domain) MSE is: 



MSE Lh = 



£ 



Appendix II 

CP-OFDM: The Class of Equispaced and Equal Pilot Tones Achieves the Global Minimum MSE 

In Appendix H] we proved that the class of equal training symbols is globally optimal for the CP-OFDM sparse 
preamble, when the pilot tones are equispaced. We now prove that this is also a globally optimal solution. 

We will rule out the possibility that nonequispaced pilot tones can yield a lower MSE than the class of equispaced 
and equal training symbols. Consider again the MSE expression (f38b . We want to minimize this, subject to the 



constraint Ylt^o 



M 



x^F mxijFmxv x Li — £■ I* i s obvious that since the CP is a wasted part of energy, 



the MSE would be minimized if somehow we were able to collect all the energy of the CP and put it in the useful 
part, namely the first sum of the constraint. To obtain the minimum MSE we can possibly imagine, we consider 
the relaxed problem^ 

MSE L „ (42) 



X 



mm 



L h* L h 

Lh-1 
S.t. Y, \ x i, 
m=0 



< £ 



(43) 



8 This is a "genie-aided" problem, i.e., a problem that is unrealistic in practice and only a genie can help us to obtain, since it would lead to 
the minimum possible achievable MSE. 
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However, it is known [19], [3] that the optimal solution for this problem is the sparse preamble of Lh equispaced 

T 2 

and equipowered pilot tones. The minimum achievable MSE is h £ . This is also achieved by a sparse preamble of 
equispaced and equal pilot tones in the previous appendix. Therefore, we only need to verify that this is the unique 
class of sparse preamble vectors that achieve this minimum MSE. This is equivalent to proving that the class of 
sparse preambles with equispaced and equal pilot tones is the only one that zeroes the CP energy. 
The submatrix of F Mxv F^ Xv involved in the evaluation ofxf h F Mxv F^ xv x Lh is given by (f 

where Xj, h is now any set of the form {io + kM/Lh\ k = Q, 1, . . . , Lh — 1} with iq = 0, 1, . . . , M/Lh — 1. This 
submatrix has a very special form: 

Lemma 2: All diagonal entries of the submatrix \F Mxv F^j xi 
its off-diagonal entries equal -1. 

Proof: The general entry of the above submatrix is given by: 



are equal to v = Lh — 1, while all 



Mxv r Mxu J . . — e h — 2-~i 

h3 1=1 1=1 

where i = + kij^- and similarly for j. Obviously, for i = j, (F mxi/^mxv)^ = — 1, i = 0, 1, . . . , Lh — 1. 
If i j, and setting k — kj — hi, 



L h -1 L h -1 

* Mxv r Mxu ) — 7, e — ?_, e 

J 1=1 1=1 



with the assumptions made previously for M, Lh- But: 

0= e J ^ M = l+(F Mxiy F M: 

1=0 

hence 

[ F Mxu F Mxu) = — 1) i^J 

■ 

The question now concerns the type of vectors x L/ that vanish the term x" h F MxiJ F Mxu x L) = \\F Mxv x L \\ 
vanish. Suppose that there is such a sparse vector with equispaced and equipowered symbols Xi m = \x\e jei ™ . Then, 
there should hold F Mxu x Lh = 0, hence F Mxv F Mxv x Lh = 0. Consider, for example, the inner product of the first 
row of F Mxu F Mxv with x Lh . Then, according to the previous lemma, there should hold \x\ {ye> ei a — Y^ n= i e j6im ) 
or ve 3<>i <' — 2~2m=i & >6im ■ Taking the modulus in both sides, we should have v = |X)m=i e° 0im \ - But this can 
only happen when all the exponentials in the sum are collinear and of the same direction in the complex plane, 
i.e., when all these exponentials are equal. 

Conclusion: Among all sparse preamble vectors, it is those with equispaced and equal pilot symbols that yield 
the globally minimum MSE for CP-OFDM. 

Remark: We can alternatively prove the statements of the last two Appendices without resorting to the Lagrange 
theory. Having defined the optimization problem (ffOt-lHTTi. we can observe that the MSE achieved by the equispaced 
and equal pilot tones when they satisfy the constraint with equality is LhO 2 /£■ Without proceeding with the proof 
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of Proposition Q] we give the genie-aided problem (I42li-(l43ll which achieves the minimum MSE for any sparse 
preamble. This minimum MSE is known to be Lh<J 2 /£ [19], [3]. Thus, we only need to verify that the class of 
equispaced and equal symbols is the unique MSE-optimal class for the sparse preamble design. This is shown as 
above. 

Appendix III 

CP-OFDM: Optimal Full Preamble Vectors With Equipowered Pilot Tones 
The MSE expression (in the time domain) for the full preamble is MSEm = Ap-tr 



where C M is the estimation noise covariance matrix, which is diagonal, with diagonal entries of the form a /\x m \ > 
m = 0, 1, . . . , M — 1. To obtain this expression, we have used the pseudo-inverse yF^j xLh F 'mxLh) F MxL h 
to translate the CFR estimates to the time-domain and the fact that \Ff IxLh F MxL^j = 0-/^)^L h - Clearly, 
here we do not face an optimal placement problem. Also, F MxL h F MxL h i s an M X M matrix with its diagonal 
elements all equal to Lh- Thus, the above MSE can be written as 

t M -! 2 

<44) 

m— 

Our problem then is to minimize this MSE subject to the constraint 2~2m=o \ x m\ 2 + jj xH F Mx „Ff Ixu x < £. 
We already know that the training vector with all equal symbols is a global minimizer for that problem. In that 
case, MSEm = j&^T = = and x H F M xu F Mxu x = °- In the following, we show that there are 

also optimal full preamble vectors with simply equipowered, not necessarily equal symbols, and demonstrate ways 
of constructing them. 

7j TT . . TT . . rj 

Obviously, if x F Mxv F Mxv x = \\F M xv x\\ = 0, then x must be spanned by the first M — v columns of the 
M x M DFT matrix T . That is, it must be of the form x = ^2f^ Q Lh cnf i7 where i = 0, . . . ,M — Lh, is the 
ith column of T. We observe that simply setting x = \J £ /M f L for any i = 0, 1, . . . , M — 1 leads to equipowered 
(but unequal) symbols that minimize the MSE (cf. d44ii). Also, due to the orthogonality of the DFT vectors and 
the fact that all of them have energy equal to M, we can see that the complex numbers oti, i = 0, 1, . . . , M — Lh 



should satisfy X)f=o Lh 



a i \ — M ' 



We now give an algorithm for constructing an infinite number of such full preamble vectors, combining at most 
two of the first M — v columns of T. Denote by F 'mx(m-v) tne corresponding M x (M — v) matrix. 



Proposition 2: We can find infinitely many (M — i^)-tuples a 



a ai ■ ■ ■ um-v-1 



T 



that satisfy 



a 



and lead to x = X^o^ 1 a ifv w i tn \ x i\ — \ x \i * = 0, 1, ... ,M — 1, in the following two cases: 
First, only one of the a's, say a m , is nonzero, leading to a scaled version of the DFT column f m , and second, 
two of the a's are nonzero, say ak,a m , with phase differences ±7r/2 and \k — m| = M/2. The second case is 
only justified if Lh < 4f . 

Proof: Clearly, the a's we look for are such that F Mx ^ M _ u ^ a = jjU, where u is any AI x 1 vector with unit 



modulus entries. For such a system of equations to be consistent, u should belong to the range space of F 



Mx(M-u)' 
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i.e., u 



Mx(M-f) 



7, for some complex vector 7 



7o 7i 



7M-I/-1 



Then a 



—7. Taking 



,1/ 



the squared norm of both sides of the last equation and using the constraint on the norm of a, we obtain 



M-v-1 

E 

i=0 



1 



Moreover, taking the first entry of u and its modulus, we can write 

M-v-1 

E 7, 



i=0 



1 



The last two equations are satisfiable in the following cases: First, if one of the 7's is unit modulus, say j m , and 
the rest of them are zero. Then, ol will have only one nonzero entry, a m , with modulus \a m \ — y^. Alternatively, 
assume that only two of the 7's are nonzero, say 7fe,7 TO , with the rest of them being zero. Then, if jk has a 



modulus 7, i.e., 7fc = ^e 30 , and j m = \/l — j 2 e J ^ e± ^ \ both equations are satisfied. In that case, there will only 



be two nonzero a's, namely 



7 fc and a r , 



7 m . Hence, for the rth entry of x, we will have 



M < m 



dike 



+ a m e 



for r = 0, and also for r 7^ if 3fJ <^e J [ 
that it implies the requirement \k — m| = 4f 



M 

"° r±f ]| = 0, i.e., if 27r(fc ," f m)r ± - = =H iu..<| ;r. This can be seen 



l|. This is obviously equal to £/M 



M 



Appendix IV 

OFDM/OQAM: Optimal Sparse and Full Preambles 
Define the vector of the nonzero SFB output samples for a sparse preamble input (cf. Section IV- Al l: 

~L h 



OQAM 



Eiezj, ai,o9i,o(°) T,iex Lh ai,o9i,o{ 1 ) 



i,Ogi,o(Lg — 1) 



Clearly, for a sparse preamble, we have M/Lh > 2. We first show the following: 

Proposition 3: If M/Lh > 2, then 1 1 s oqam II 2 = EieZz, a ?0' Le -> tne ener gy transmitted for training is equal to 
the energy of the training vector at the AFB output of the associated ideal (channel- and noise-free) OFDM/OQAM 
system. 
Proof: 

L a -1 



1=0 
L„-l 



Obviously, 



and 



Lg-l 




2 


L h -1 L h -1 


s Lh II 2 - V 
6 oqamII — / / 


E a i;09i,o( l ) 




E ai ™>° E a 


1=0 






m=0 k=0 


L h -1 


Lg-l 




L h -1 


= E 


3 L,o E \si m ,o 


(0I 2 


+ E ai ^' 


m=0 


1=0 








Lg-l 




Lg-l 




E i^,o(oi 2 = 
1=0 


= E 5(0 2 = 1 

1=0 



L a -1 



E Si m ,0(0^,0(0 = °> m,k = 0,1,..., L h -1, m ^ k 



1=0 
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due to our assumption that M/Lh > 2 and a good frequency localization of g. Therefore, || Sqqam II 2 = J2m=o a f m 0> 
i m € T Lh . ■ 
The (time domain) MSE expression for the sparse preamble with Lh nonzero pilot tones in the OFDM/OQAM 



system is the same as in the CP-OFDM system, i.e., MSE^ = tr 
problem can therefore be stated as follows: 



C T , Fr.»T..F 



H 



Our optimization 



min MSE Lh (45) 

a im ,i m £l Lh 
L h -1 

s.t. ^ a\ m < £ (46) 

m— 

where we have suppressed the temporal index 0. But the solution to this problem is known. It is the class of 
equipowered and equispaced pilot tones [19], [3]. 

For the full preamble, it is easy to show that the transmit training energy is not equal to the energy of the training 
vector at the AFB output. Using our assumption on the time-frequency localization of the prototype function, the 

training energy constraint can be written as: 

M-l 

^ ] (l^ml + P x m x m —l P x m x m +l) — & ■> 

where x m — a m fie 3Vm -° , and therefore the optimization problem can be stated as: 

min MSE M = — V ^ (47) 

x m ,OT=0,l,...,M-l M ^ ' \x m + pxm-i + px rn+ i\ z 

M-l 

S.t. (l X ™| 2 + P x ^ x * m -l + famX* m +l) < £ (48) 

m=Q 

Using similar steps as in the proof of Proposition Q] we can easily show that the full preamble vector with all equal 
symbols is a minimizer of (07), d48l i and the minimum achievable MSE is Ma 2 / [£(1 + 2/3) 2 ]. Furthermore, we 
can show that this preamble is a global minimizer of the last optimization problem, but with a constraint on the 
training energy at the SFB input, i.e., with a constraint of the form: 



]T \x m \ 2 <E 



Proposition 4: A global minimizer of the problem: 

I A/-1 2 

min MSE M = — V -, ^ ^ pr (49) 

x m , m =o,i,...,M-i M ^ ' \x m + px m -i + /3x m +ir 

m—0 1 'ii 

M-l 

s.t. i x ™i 2 ^ £ ( 50 > 

m=0 

is the full preamble with equal symbols. 

Proof: To prove the statement of this proposition we can initially show that the equal symbols is a local 
minimizer of our optimization problem via Lagrange theory and then verify that the equal symbols lead to the 
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minimum possible MSE. However, if we show that the equal symbols achieve the lowest MSE, the step associated 
with the Lagrange theory in unnecessary!^. 

We can divide the class of symbols into the following subclasses: 

la) Equal symbols, i.e., symbols of the same magnitude and phasJ"! 

lb) Symbols of equal modulus and different phases (at least one symbol with different phase from the rest of 
the symbols in the preamble vector). 

2a) Symbols of different modulus but of the same phase (at least one symbol with different modulus than the 
rest of the symbols in the preamble vector). 

2b) Symbols of different modulus and phase (at least a symbol different from the others in the preamble vector). 

The subclasses la), lb) subdivide the general class of equipowered symbols and subclasses 2a), 2b) the general 
class of nonequipowered symbols. First, we will show that between subclasses la) and lb), la) leads to an equal 
or lower MSE than that of lb). The same holds for the subclass 2a), when it is compared with 2b). 

Comparison of la), lb) subclasses: We consider a vector of equal symbols xq = x\ = ■ ■ ■ = Xm-i = x = \x\e 3 ^ 
and an arbitrary vector of the form x = \x\e 3 ^°, X\ — |x|e J<?il , . . . , xm-i — \x\e 3< ' >u - 1 . For the case of equal 

2 2 

symbols, the arbitrary term ^ +f3x — a 1+ p x |2 takes the value \ x \ 2(i +2 0)2 ar, d if we assume that the constraint is 
satisfied with the equality, then this value becomes g7^§^p ■ In the case of unequal symbols, the modulus is again 
£/M since the constraint is phase blind, thus the maximum value that can be taken by any such term is gn^^p ■ It 

2 

is obvious that in any way the phases Sq, d>±, . . . , 4>m-i are chosen, there is at least one term , — -r-r, — - „ r? 

for the unequal symbols that is greater than or equal to g7^§^p (essentially, we can not achieve in all cases triplets 
of numbers x m ^i,x m ,x m +i that are collinear and of the same directionality in the complex plane). Therefore: 

MSE lQ < MSE 1& 

In the same way: 

MSE 2a < MSE 26 

To finish this proof, we have to show that MSE la < MSE 2a . For the subclass la) we can write: 

1 x-^ a Ma 



MSE lQ = — V 
M ^ 



M \ Xm + /3:Em " 1 + ^™+il 2 YZZI km + I3x m ^ + /3x m+1 \ 2 

Ma 2 
£(l + 2/3) 2 

using the Arithmetic-Geometric-Harmonic (AGH) mean inequality. For the subclass 2a), we have: 

Ma 2 



\ISE 2a = — V r 



9 Note that this holds also for the proof of Proposition 1. I.e., we can alternatively combine the results in the first two appendices, discarding 
the Lagrange theory step and simply verifying that the equal and equispaced symbols achieve the lowest MSE, which is equal to L^a 2 /£. 
l0 We refer to the interval [0, 2tt) for the phases, since for a phase 4>o 6 [0, 27r], the phase 0o i 2k-K, k £ Z leads to the same phasor e J ^°. 
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If we show that: 

M-l 



m=0 

or 

M-l 



]T (Cm + /3Cm-l + Km+lf < E{1 + 2/3) 2 



m=0 



we will be done. Here, Cm — \x m \ and the last equation is obtained due to the equal phases of the symbols. 
For the last expression, we have: 



E™=0 (Cm + IKn-1 + /3Wl) 2 = £ C +/? 2 E Cm-1 + ^ £ C 



M-l M-l M-l 

2 

m+1 

ra— m— m— 



M-l M-l M-l 

+ 2/3 £ CmCm-1 + 2/3 E CmCm+1 + 2/? 2 £ Cm-lCm+1 

m— m— m— 

We can consider that the index m — 1 for m = equals M — 1 and correspondingly, the index m + 1 for m = M — 1 
equals due to the periodicity of the discrete time spectrum. Therefore: 

M-l M-l M-l 

Cm = Cm-1 = Cm+1 = & 

m—0 m—0 m—0 

For the same reason: 

M-l M-l 

^ ^ CtnCm — 1 ^ ^ CmCm+1 
m—0 m—0 

We set Cm-i = Cm in the last equation. We place the numbers Cm and £ TO to the main diagonals of two diagonal 
matrices Z and H, respectively. The following lemma [8, p. 183] will prove useful: 

Lemma 3: Consider any two matrices A,B e c m xn and their singular values, di(A), <7i(B), i — 1,2, ... ,q, 
with q = min{m, n}, arranged in a descending order. Then: 

(a) The following inequality holds: \tv(A H B)\ < J^Ui ^i{ A )^i{ B )- 

(b) There exist unitary matrices P 1 and P 2 such that max{|tr(P 1 A H P 2 S)| : P 1 e C nxn ,P 2 € C mxm arc unitary} = 
Using the last lemma, we write: 

M-l M-l 

tr(ZS)= ECmCm-l< ECm=f 
m—0 m—0 

Using the same trick, we can show that: 

M-l M-l 

E Cm+1 Cm-1 < E Cm = £ 
m—0 m—0 

Therefore, 

M-l 

E(Cm + /3Cm-i+/3Cm+i) 2 < £ + 0> E + p 2 £ + A0E + 2P 2 £ 

m=0 

= £(l + 2(3) 2 

and our proposition is proved. ■ 
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Appendix V 

OFDM/OQAM: MSE with Postprocessing of the Full-Preamble-Based Estimates 

To derive the MSE expression, we first have to determine the covariance of the noise at the AFB output. The 
mth noise component at the output of the AFB of the OFDM/OQAM system corresponding to the preamble vector 
is given by 7] m Q — Y^f=o ^(Offm o(0> wnere W (J) i s tne noise at the receiver front-end, assumed additive and 
white with zero mean and variance a 2 . We then have: 

E (Vm,oVk,o) =J2J2 E [ w ( l ) w *i r )] 9*n,o( l )9kA r ) = 0-2 <VSm,o(0sM( r ) = a2 <4,o(0.9fc,o(0 

I r I r I 

For m = k, this is equal to a 2 . According to our assumptions on the good localization of the pulse in frequency, 
and using the fact that the phases, ip m .o, of the symbols are equal, we can see that the cross-correlation term above 
takes the value a 2 (3 for I = m — 1, m+ 1 and is zero elsewhere. Therefore, the noise covariance matrix at the AFB 
output can be expressed asF*l 

1 p ••• ±p 

, P i P ■■■ o 

C M = a 2 



±P 



ft 



B 

When computing the initial CFR estimates as in ( flOl ), the noise covariance matrix becomes C' M = £ oiiAM^ 1+2 ^ 2 B, 
and hence the corresponding MSE is Ma 2 / [£° QAM (1 + 2/3) 2 ] as in ( fLST l. After the processing described in Sec- 
tion IV-CI we come up with the final CFR estimates, with MSE given by: 



MSE 



£2 QAM (l + 2ft) 



-tr 



F]MxL h (^MxLh^MxLh) ^MxL h BF MxLh (^Ff IxLh F 



MxL h 



MxL h 



£OQAM (1 + 2/ j )2M 



— tr (BF MxLh F I Mxl , 



? MxL h Ffl: 



L has all its main diagonal entries equal to L^. For large M, it is easy to see that its entries 
immediately above and below its main diagonal can also be well approximated by Lh- This leads to: 

MSE PS r>r>AAj 

£° QAM (l + 2/3) 
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